Resolving the Sub-AU-Scale Gas and Dust Distribution in FU 

Orionis Sources 

J.A. Eisner 1 ' 2 & L. A. Hillenbrand 3 



ABSTRACT 



We present Keck Interferometer observations of the three prototypical FU 
Orionis stars, FU Ori, V1057 Cyg, and V1515 Cyg. With a spatial resolution of 
a few milli-arcseconds and a spectral resolution of A/AA = 2000, our near-infrared 
observations spatially resolve gas and dust emission extending from stellocentric 
radii of ~ 0.05 AU to several AU. We fit these data with accretion disk models 
where each stellocentric radius of the disk is represented by a supergiant-type 
stellar emission spectrum at the disk temperature. A disk model is consistent 
with the data for FU Ori, although we require some local asymmetry in the disk. 
For V1057 Cyg the disk model does not fit our data well, especially compared 
to the fit quality achieved for FU Ori. We speculate that a disk wind may be 
contributing substantially to the observed near-IR emission in this source. The 
data for V1515 Cyg are noisier than the data obtained for the other two objects, 
and do not strongly constrain the validity of an accretion disk model. 

Subject headings: stars:pre-main sequence — starsxircumstellar matter — stars:individual(FU 
Ori, V1057 Cyg, V1515Cyg) — techniques: spectroscopic — techniques: interferometric 



1. Introduction 



The class of FU Orionis objects was initially defined by s udden, large incre ases in 
optical-wavelength fluxes in FU Ori, V1057 Cyg, and V1515 Cyg (jHerbiglll966l . Il977l ). After 
these brightening events the optical spectra of these objects resemble those of low surface 
gravity (i.e., supergiant) stars. However for V1057 Cyg, pre-outburst spectra show that th e 
progenitor resembled a higher surface gravity, typical T Tauri type star (e.g.. iHerbigi Il958h . 



1 Steward Observatory, University of Arizona, Tucson, AZ 85721 
2 Alfred P. Sloan Research Fellow 

3 Astrophysics Department, California Institute of Technology, Pasadena, CA 91125 



-2 - 



While there are n ow many more members of the FU Ori class (e.g., iReipurth fc Aspinl 12010 
Miller et al.ll201ll ). the three prototypes remain very well-studied. 



FU Ori outbursts are believed t o result from periods of enhanced accretion in circum 



stellar disks arou nd young stars (e.g.. lHartmann fc Kenyonlll985t lHartmann fc Kenvon 



Zhu et al.l 120071 ). Outbu rsts may be triggered b y thermal instabilities (IBell &: Lin 1994), 



gravitational instabili ties ( iVorobyov fc Basull2005l ). a combination of the two ( lArmitage et al. 



1996 



2001 



Zhu et al.ll2009l). or by interactions of a disk with a planet or nearby stellar com panion 



(ILodato & Clarke! 12004 IClarke et al.l 120051 : IReipurth fc Aspinl 12004 |Pfalzner!l2008f l. Spec- 
troscopic observations of FU Ori stars have al so been explained in the context of rapidly 



rotating stars near the edge of stability (e.g., iHerbig et al.l 120031 ) . However spectroscopic 
data at a range of wavelengths exhibit various features that are more consistent with disks 
than stellar photospheres, including highly broadened lines and poor c ross-correlation sig- 
nificance with single-temperature, fast-rotating giant star spectra (e.g., iGreene et al.ll2008l . 
Hillenbrand et al. in prep.). 

The disk model is strongly supported by interferometric observations that resolve the 
circumstellar emission at near-IR and mid-IR wavelengths on scales > 100 R . Near- and 
mid-IR observations of FU Ori itself can be explained in the context of a geometrically thin 
accretion disk extending t o stellocentric radii as small as 5 . 5 R , inclined 55° with respect 
to the plane of the sky (IMalbet et al.l 120051 ; lOuanz et al.l 120061 ) . The high inferred disk 



temperature implies an accretion rate > 10 5 M yr 1 , substantially higher than for typical 
T Tauri stars. The disk model for FU Ori supported by the interferometric d ata is also 



comp atible with spectrophotometric data across a broad wavelength range (e.g.. IZhu et al. 
2008h . 



However the disk model may not be the whole story, and other components like extended 
envelopes or outflows may be important in explaining the data for other FU Ori sources. 
The sizes of near-IR emitting regions in V1057 Cyg and V1515 Cyg are larger than expected 
from pure disk models, su ggesting the presence of scattered light from extended structures 
(IMillan-Gabet et al.l 120061 ) . Indeed, these objects also display e xcess long- wavelength emis- 



sion that can not be reproduced w ith geometrically thin disks (IKenyon fc Hartmannl 11991 



Green et al. 


2006: 


Zhu et al. 


2008) 



To explain these data, models with both disks and 



circumstellar envelopes have been constructed. 

A key prediction of accretion disk models is differential rotation. A Keplerian disk 
rotates faster at smaller stellocentric radii, and shorter-wavelength spectral lines are expected 
to have broader linewidths. Indeed, such behavior has been observed in a number of FU 



Ori objects (e.g., lHartmann et al.l 12004 : IGreene et al.l 120081 : IZhu et all 120091 ) . While often 



difficult to reconcile with available interferometric data, alternative explanations for these 
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line profiles have been proposed (e.g.. iHerbig et al.ll2003l ). Spatially resolved near-IR spectra 
can determine the location of the gas directly, and provide a direct test of the rotating disk 
model. 

Recent enhancements to the Keck Interferometer (KI) now enable spatially resolved 



spect roscopy of the near-IR emission from circumstellar disks ( lEisner et al.ll2010bl ; iPott et al. 



2010l ). Here we present KI observations of the three prototypical FU Ori objects: FU Ori, 
V1057 Cyg, and V1515 Cyg. Our data probe the entire K-baiad and are sensitive to emission 
from the Br7 transition of hydrogen and the ro- vibrational overtone transitions of CO. We 
determine the relative spatial distributions of line and continuum emission in these systems, 
providing critical tests of circumstellar disk models. 



Observations and Data Reduction 



KI is a fringe-tracking long b aseline near-IR Michelson interferometer combinin g light 
from the two 10-m Keck apertures (jColavita fc Wizinowich!l2003l ; IColavita et al.ll2003l ). Each 
of the 10-m apertures is equipped with a natural guide star adaptive optics (NGS-AO) system 
that corrects phase errors caused by atmospheric turbulence across each telescope pupil, and 
thereby maintains spatial coherence of the light from the source across each aperture. Optical 
beam-trains transport the light from each Keck aperture down into a tunnel connecting the 
two Kecks and to a set of beam combination optics. 

We used the "self phase referencing" (SPR) mode of KI, implemented as part of the 
ASTrometric and phase- Referenced Astronomy (ASTRA) program (see Woillez et al. 2011 
for a description of the SPR mode). SPR allows deep integrations in a spectroscopic channel, 
which enables higher dispersion without loss of signal-to- noise. We use a grism with A/ A A = 
2000 to disperse the light. Note, however, that spectra are fully Nyquist sampled only at a 
resolution of ~ 1000. Near-IR spectral lines in FU Ori sources may be as broad as ~ 100 



km s 1 (e.g.. iGreene et al.ll2008 



Biscaya et al.lll997l ). and thus the spectral resolution of KI 



is reasonably matched to the observed systems. 

Neon lamp spectra and Fourier Transform Spectroscopy are used to determine the wave- 
length scale for each night of observed data. While the entire i^-band falls on the detector, 
vignetting in the camera leads to lower throughput toward the band edges. The effective 
bandpass of our observations is approximately 2.15 to 2.36 jum. 

We obtained Keck Interferometer (KI) observations of V1057 Cyg on 2009 July 15, of 
V1515 Cyg on 2010 July 21, and of FU Ori on 2010 November 25 (all dates UT). On each 
night we also observed spatially unresolved stars of known spectral types, used to calibrate 
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the instrumental response. The calibrator stars are HD 192985 and HD 195050 for V1057 
Cyg; HD 194206 and HIP 102667 for V1515 Cyg; and HD 37147 and HD 42807 for FU Ori. 

We measured fluxes, squared visibilities (V 2 ), and differential phases (A</>) for our targets 
and calibrator stars in each of the 330 spectral channels across th e if -band provided b y the 
grism (Figures [TH3]). The data calibration is described in detail in lEisner et al.l ( 120 1 Obi ) . but 
we provide a basic sketch here. 

Normalized fluxes are calculated by dividing observed source and calibrator spectra and 
then multiplying b y an appropriate stella r template for the calibrator (given by a Nextgen 
model atmosphere; lHauschildt et al.lll999l ). The absolute flux level and slope of the spectra 
are then matched to near-IR photometry from 2MASS. Error bars for the calibrated fluxes 
reflect channel-to-channel uncertainties, but do not include uncertainties in the overall nor- 
malization. 

Squared visibilities for targets are divided by V 2 measured for calibrators. As for the 
fluxes, error bars reflect only channel-to-channel uncertainties. We do not include the un- 
certainty in the overall normalization of the V 2 . 

Differential phases-i.e., the phase in each spectral channel relative to some fiducial 
channel-are calibrated by subtracting the A<p measured for the calibrator stars. The broad- 
band slope and curvature of the measured differential phases are altered by instrumental 
effects. The slope is altered by errors in group delay tracking, and curvature in the phase 
versus wavelength may be introduced by atmospher ic water vapor tu rbulence and dispersion 
and refraction during internal beam-transport (see iPott et al.ll2010l ). Hence we can not at- 
tribute any broadband slope or curvature in the data to the observed source. We therefore 
remove the slo pe and curvature seen across our bandpass from the differential phase data, 
as described in IPott et al.l (120101 ) . 



3. General Features of the Data 

All objects in our sample show strong absorption in the spectral regions corresponding 
to the CO ro- vibrational bandheads, v = 2 — >• 0, 3 — >■ 1, and 4 — >■ 2 (Figures HHS]). The rest 
wavelengths of these three bandheads are 2.2936, 2.3227, and 2.3527 /mi, respectively. 

FU Ori and V1057 Cyg also show dips in V 2 at the wavelengths corresponding to these 
spectral features. V1515 Cyg may show similar dips relative to the V 2 seen in adjacent 
channels corresponding to continuum emission, but these dips are not significant given the 
noise in these data. 
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Lower V 2 generally means more resolved emission. However, because the CO appears 
in absorption, the dips in V 2 actually indicate that the CO distribution is more compact 
than the continuum emission. This can be understood by considering a compact, hot region 
of the disk that includes CO absorption, and a cooler, continuum-emitting region at larger 
radii. At wavelengths where CO absorbs emission from the hot, compact component, the 
flux-weighted size of the two regions moves to larger stellocentric radii. Thus, the compact 
CO absorption leads to a larger measured size at these wavelengths, and hence a dip in V 2 . 

The observed differential phases indicate a centroid offset between the continuum emis- 
sion and the CO absorbing material. This offset is most obvious for FU Ori, and marginally 
significant for V1057 Cyg. For V1515 Cyg the data are too noisy to argue the presence or 
absence of such an offset. Differential phase offsets in the CO bandheads indicate that the 
annuli of the disk that have appropriate temperatures for CO absorption in their spectra 
must have centroids that differ from annuli that are too hot for CO to exist or cool enough 
that they produce continuum emission. 

The spectra for all objects show an absorption feature at the central wavelength of Br7, 
2.f662 /im. In addition, the spectra show hints of some emission on the red side of the line. 
These spectral sha pes are consistent with P Cygni profiles seen in other hydrogen lines (e.g., 



Herbig et al.ll2003h 



One object-V1057 Cyg-also shows a significant dip in V 2 in the Br7 spectral region. 
Since the line is predominantly seen in absorption, the same argument applied to CO above 
suggests that the Br7 absorption occurs in more compact regions than the continuum emis- 
sion. The Br7 feature is more prominent in the spectrum of V1057 Cyg than in FU Ori or 
V1515 Cyg; the higher line-to-continuum ratio of Br7 in this object may explain why the 
V 2 dip is only seen here. 

The CO bandhead features are much broader in VI057 Cyg than in other objects in 
our sample. These broad features do not appear to be due to instrumental effects. During 
the night of our V1057 Cyg observations, we observed known CO-emitting objects immedi- 
ately before and af ter V1057 Cyg. Neither of these objects show such broad features (see 



Eisner et al.ll2010a[ ). We also verified that the appearance of the CO features in V1057 Cyg 
persists across multiple observations and for different calibrator stars. Thus the broadening 
seen in the CO lines of V1057 Cyg appears to be intrinsic to the source. 
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4. Modeling 



In this section we attempt to create a physically-motivated cirumstellar disk model to 
explain the measured fluxes, squared visibilities, and differential phases. In our modeling 
we assume that the emission from our targets arises entirely from circumstellar material. In 
previous modeling of FU Ori sources spectrophotome tric data has be en fitted well with no 
contribution from a compact stellar photosphere (e.g.. lZhu et al.ll2008l ). We therefore do not 
require any correction for unresolved flux from the central star. 

We begin with a geometrically thin disk model. We assume the disk spans a range of 
radii from R in to R out , and has a temperature profile T = T in (R/R in ) a . We divide the disk 
up into a series of annuli and assume the emergent spectrum of each annulus is that of a 
supergiant star. For FU Ori stars this is a reasonable assumption since they are heated 
primarily from accretion energy released near the midplane. Disk atmosphere models have 
been constructed specifically for FU Ori systems, and these do provide bett er fits to detailed 
spect rophotometric data than models based on stellar atmospheres (e.g., lHartmann et al. 
2004 ). We are, however, not concerned with this level of detail here, and we deem our 
models based on giant star atmospheres to be sufficient. 

For annula r temperatures between 1700 and 10000 K, we use emergent spectra for stars 
with \ogg = (jHauschildt et al.l 119991 1 . (In practice, the maximum temperature is typically 
between 4000-6000 K.) For T < 1700 K, we assume that the emission is a simple blackbody. 
The emergent spectrum is rotationally broadened by the Keplerian velocity expected for each 
annulus, v sin i = y/GM/Rsini. We take M = 0.5 M for our obje cts; this is compa tible 
with the dynamical mass of 0.3 ±0.1 M for FU Ori determined by IZhu et al.l (120091 ) and 
the pre-outburst spectral type of V1057 Cyg. The spectrum of the disk model is given by 
the sum of fluxes over all annuli. 



The visibility of each annulus is given by the visibility of a uniform ring (e.g. JEisner et al. 



20041 ) . The visibilities for the entire disk are computed as the flux- weighted average of the 
annular visibilities: 

^annulus ^annulus \ 
y j -^annulus / 



V. 



disk 



(1) 



The differential phase is related to the centroid offset. For each annulus, the differential 
phase is 

A0 annu l us ~ 27Tr uv A$ annu i us , (2) 



where r uv is the uv radius (see, e.g., Eisner et al.l 12004 ) and A^^ius is the centroid offset. 
This approximation is valid as long as the emission is not well-resolved at the angular res- 
olution of the observations. The differential phase signal for the entire model is (see, e.g., 
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Eisner et al.ll2010bf ) 



• , y j ^'/'annulus^'annulus-^annulus / \ 

A 0disk = f^tt p • (3) 

/ / ' annulus^annulus 

As described above, the slope and curvature of the observed, broadband differential phase 
are altered by instrumental effects. We therefore remove broadband slope and curvature 
from both observed and modeled A0. 

The simple, geometrically thin disk model considered above would not produce centroid 
offsets. However if we allowed the disk to be somewhat asymmetric, certain disk radii could 
have centroids offset from those of other disk radii. We allow for the possibility of such 
asymmetry by including a centroid offset that changes linearly with disk radius. This aspect 
of the model is described by a single free parameter that gives the amplitude of this offset. 
We denote this amplitude by A# max . 

To account for the possibility that some of the near-IR emission arises from more spa- 
tially extended regions-f or example in remnant envelopes (as implied by the results of 



Millan-Gabet et al.l l2006l )-we include a scattered light contribution in our models. While 
there may also be thermal emission from a more extended component, the cooler tempera- 
tures at large stellocentric radii would produce a very small fractional flux in the near-IR. 
We assume that the extended, scattered emission has the same spectrum as the disk, and is 
sufficiently extended so that is is completely resolved (V 2 = 0) in our observations. Writing 
the fractional scattered flux, / = -F scat /Fdi S k, the model fluxes and visibilities are modified 
to 

-Fmodel = -Fdisk(l + /), (4) 

VLhl = + f)~ 2 - (5) 

Finally, we include two scaling factors in the model that take into account additional 
uncertainties not reflected in the error bars for the data. The first is a constant scaling in the 
overall flux level, Xfl ux . This scaling accounts for uncertainties in the flux calibration beyond 
the assumed channel-to-channel errors, including uncertainties in the 2MASS photometry 
used to normalize the spectra, or source dimming between the photometric and interfero- 
metric observations. We also include a constant scaling in the overall V 2 normalization, X v i , 
to account for overall calibration uncertainties beyond the channel-to-channel uncertainties 
reflected in the plotted error bars. These two scale factors may also absorb uncertainty in 
the assumed source distances (Tabled]). It turns out that these scaling factors are near unity 
in our best-fit models, and hence are not crucial to the fitting process; we include them for 
completeness. 

For FU Ori, the disk inclination and position angle have been measured using previous 
interferometric observations. In the near-IR the inclination and position angle were estimated 
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as ~ 55° and 47°, respectively ( IMalbet et al.l 120051). In the mid -IR, the inclination and 
position angle were found to be ~ 55° and 100° (jOuanz et al.l 120061 ). We adopt the near-IR- 
determined values here. For V1057 Cyg and V1515 Cyg, no such geometrical measurements 
are available, and so we leave inclination and position angle as free parameters in our models. 
We can not, however, constrain these well with the limited uv coverage of our observations. 

Given the large numbers of free parameters, we fit our models to the data using a 
non-linear gradient-fitting algorithm based on the Levenberg-Marquardt method. Such al- 
gorithms may fall into local minima in the x 2 surface, and so initial parameter guesses and 
step sizes for gradient searches influence the results. We have selected step sizes that ensure 
that all free parameters are varied during the fitting process. However we are hesitant to 
claim that the "best-fit" models we determine are the true, global best fits. 

To improve the signal-to-noise ratio in the fitting, we only use data with wavelength 
between 2.15 and 2.36 /im to constrain the models. Data at the band edges are substantially 
noisier because of instrumental vignetting and contributions from ambient water vapor (Pott 
et al. 2010; Woillez et al. 2011). 



5. Modeling Results 
5.1. FU Ori 

Our best-fit disk model for FU Ori is shown in Figured! The model fits quite well in this 
wavelength region (reduced x 2 = 0.5). While the modeled differential phases do not match 
the observations perfectly, this is likely due to imperfectly subtracted slope and curvature 
(which are instrumental effects). Indeed, residual curvature of the differential phase data 
forces the best-fit model toward a flat differential phase versus wavelength. 

If we consider the differential phase data only over the narrower spectral region contain- 
ing the CO bandheads (where broadband slope and curvature effects are less important) we 
can easily construct a model that agrees with the data well. In particular, we use the best-fit 
model parameters for FU Ori listed in Table [2J increase the value of fl max , and subtract a 
constant phase offset. Note that # max has no impact on the fluxes or V 2 values and so we are 
free to modify its value in order to better-fit the differential phases. The dashed line in Fig- 
ure [1] shows that a model with a large centroid offset, # max = 0.8 mas, produces a differential 
phase signature that agrees well with the observations in the CO bandhead region. 

The differential phase signature seen for FU Ori thus can not be explained with a 
simple, geometrically thin disk model. Rather we require a centroid offset that depends on 
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stellocentric radius. A disk warp can obscure part of the outer regions of a disk while leaving 
the inner regions unobscured (or vice versa). The centroid offsets implied by our data suggest 
that the FU Ori inner disk may be significantly warped. In our model this warp occurs over 
stellocentric radii spanning 0.07 AU to 0.49 AU. 

Disk winds may also cause centroid offsets that vary as a function of stellocentric radius 



^e.g., iCalvet et al.lll993l ; lHartmann fc Calvet| 1995). Since all o f the sources in our sample 
are known to be losing mass via winds (e.g.. lCroswell et al.lll987l ). this may be an alternative 
explanation for the differential phase signatures seen in the data. In any event, the large 
differential phases observed across the spectral region including CO bandhead absorption 
indicate substantial geometric asymmetry in the underlying emission from FU Ori. 

The outer radius derive d for FU Ori is ~ 0.5 AU, similar to the ~ 0.5-1 AU value derived 



in the the SED modeling of IZhu et al.l ( 120071 ) . This likely rep resents the outer radius of the 
hot, high-accretion rate, inner disk component. As argued in IZhu et al.l ( 120071 ). this implies 
that the high accretion rate in FU Ori is probably not caused by pure thermal instability, 
but rather gravitational or MRI instability. Note that a larger R out value is not required to 
fit the broadband SED, and our model fits photometry out to A ~ 10 well. To fit longer 
wavelength pho tometry requires an additional model component such as a flared disk (e.g., 
Zhu et alJliooi ). 



5.2. V1057 Cyg 



For V1057 Cyg the disk model does not fit the CO bandhead features well (Figure [2]). 
The reduced x 2 residuals of the fit to our if-band data are 1.9. The discrepancy is due to the 
extremely broad CO overtone absorption features. For an assumed 0.5 M Q central object, 
the Keplerian velocities do not exceed ~ 100 km s" 1 , and hence the modeled CO bandheads 
are not broadened significantly. 

A higher inclination would yield higher v sin i values, but would also lower the modeled 
flux and make it difficult to fit the observations. Moreover even for i = 90° the synthetic 
spectra are still not as broad as the data. A higher central mass can also increase the 
Keplerian velocities. However M* > 50 M Q is required to produce the ~ 500 km s _1 
broadening seen in the data. A stellar mass this high is n ot (remotely) compatible with the 
pre-outburst, K-spectral type determined for V1057 Cyg (]Herbiglll958h . 



Such high velocities could, however, be created by a rotating dis k wind (e.g..lCalvet et al. 



1993). V1057 Cyg is known to possess a variable, high- velocity wind (lHerbigll2009t iCroswell et al 
19871 ). If this wind dominated the near-IR absorption during the epoch of our observations, 
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it may explain the broad CO absorption features in our data. While prev ious observations 



of VI 05 7 Cyg found the C O bandheads to be narrower than inferred here (ICarr et al.lll987 



Hartmann &: Kenvonlll987n . the CO overtone absorption spectra show substantial variability 
( IBiscaya et al.lll997l ). When the wind is more dominant, the accretion luminosity generated 
in the disk midplane may be absorbed by CO (and other matter) with radial velocities 
substantially higher than Keplerian. 

Rotating disk-wind models are substantially more complex than the simple disk models 
considered here. Rather than creating an accurate wind model, we demonstrate the overall 
features of such a model with a simple modification of our baseline disk model. Specifically, 
we scale up the Keplerian rotation profile to account for the possibility of higher radial 
velocities due to outflowing motions. 

Figure |2] shows that a model with substantially higher radial velocities does indeed fit 
the data better than the simple disk model. Here we scaled the disk velocity structure up 
(arbitrarily) by a factor of 30, so that the CO lines are broadened by vsini ~ 500 km 
s _1 . This model fits the data substantially better than the Keplerian disk model considered 
above; the high-velocity model produces a reduced x 2 value of 1.3. While the fit deviates 
from the data for the CO v = 4 — >• 2 transition, this deviation may be due, in part, to lower 
signal-to-noise of the data toward the band edge. 

Our best-fit disk model also fails to reproduce the V 2 dip in the Br7 spectral region 
observed for V1057 Cyg. Because the inner disk temperature is ~ 4000 K for our best-fit 
disk model, there is essentially no Br7 absorption in the model spectra. To explain the 
Br7 feature seen in the data would require an additional model component with a higher 
temperature. This component would also have to be reside at very small stellocentric radii 
(< 0.1 AU, which is more compact than the distribution of CO-absorbing material) to 
produce the observed dip in V 2 . Across the K-band, a simple disk model is not an obvious 
match to our observations of V1057 Cyg. 



5.3. V1515 Cyg 

The V 2 and A0 observed for V1515 Cyg are noisy compared to the data for the other 
two objects in our sample. V1515 Cyg is more than a magnitude fainter than V1057 Cyg, and 
~ 3 mags fainter than FU Ori, and so more noise is expected. V1515 Cyg is near (but within) 
the brightness limit of the SPR mode, beyond which signal-to-noise degrades significantly 
(Woillez et al. 2011). However, lower throughput near the band edges means that the signal- 
to-noise is degraded more in the spectral region including CO. This low signal-to- noise gives 
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rise to the un-physically large V 2 (i.e., normalized V 2 > 1) seen for A > 2.33 /xm. 

The fluxes for V1515 Cyg appear reasonably consistent with our simple disk model 
(Figure [3]). The V 2 and A(f> are not very consistent with the data, as seen in the high reduced 
X 2 value for the fit (Table [2]). The large \ 2 value is due, at least in part, to underestimated 
uncertainties. The large V 2 seen at the long-wavelength edge of the band are unphysical, 
and so the plotted error bars are clearly underestimated. Over most of the band the model 
appears compatible with the noisy V 2 and A0 data. 

6. Conclusions 

We presented near-infrared, spatially resolved, spectroscopic observations of FU Ori, 
V1057 Cyg, and V1515 Cyg, obtained with the Keck Interferometer. With R = 2000 across 
the .fT-band, we spectrally resolved individual CO bandheads as well as the Br7 transition of 
hydrogen. Our measured fluxes, squared visibilities, and differential phases were compared 
with simple accretion disk models. 

A disk model is consistent with the data for FU Ori, although we require some asym- 
metry in the inner ~ 1 AU of the disk to explain the measured, non-zero, differential phases. 
We speculate that the inner disk of FU Ori may be warped. The disk model also appears 
reasonably compatible with the data for V1515 Cyg, although these data are noisy and so 
the constraints are weak. 

For V1057 Cyg the disk model does not fit our data well. The observed CO bandhead 
features are broader than expected for a disk in Keplerian rotation. The disk model also 
fails to reproduce compact Br7 emission seen in the data. We speculate that a rotating, 
high-velocity, inner disk wind may be contributing substantially to the observed near-IR 
emission in this source, at least at the epoch of our observations. 

Overall this study suggests a more complex explanation for the near-IR observations 
of three canonical FU Ori stars than the unified picture of a simple accretion disk model. 
In addition to a rotating disk undergoing viscous accretion, outflow is also likely to be an 
important element in the interpretation of radiative and kinematic data. While the namesake 
object, FU Ori, is consistent with a simple accretion disk model, V1057 Cyg is not. It is 
possible that V1057 Cyg is undergoing an epoch of enhanced outflow, and that other FU 
Ori objects experience similar episodes. However further observations, at additional epochs, 
are needed to test this hypothesis. 

Data presented herein were obtained at the W. M. Keck Observatory, from telescope 
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Fig. 1. — Fluxes, squared visibilities, and differential phases observed for FU Ori (data 
and la uncertainties shown with gray shaded regions). We also plot the fluxes, V 2 , and 
A<j) predicted by our best-fit disk model (solid curves). The best-fit model predicts a flat 
differential phase across the observed bandpass, and misses the significant deviations in A(j) 
seen in and out of the CO bandheads. A model exhibiting a centroid offset with A# max = 0.8 
mas (dashed curve) fits the differential phases in the CO bandhead region well. 
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Fig. 2. — Fluxes, squared visibilities, and differential phases observed for V1057 Cyg (data 
and la uncertainties shown with gray shaded regions). We plot the fluxes, V 2 , and 
predicted by our best-fit disk model with solid curves. We also show the predictions of a 
model where we scaled the velocities up so that the CO lines are rotationally broadened by 
~ 500 km s^ 1 (dotted curves) to illustrate how a disk wind might appear different than a 
rotating disk. The model with higher velocities produces a reduced x 2 °f 1-3, compared to 
f .9 for the baseline disk model. 
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Fig. 3. — Fluxes, squared visibilities, and differential phases observed for V1515 Cyg (data 
and la uncertainties shown with gray shaded regions). We also plot the fluxes, V 2 , and A0 
predicted by our best-fit disk model (solid curves). Note that the V 2 values greater than 
unity, seen for A > 2.33 /im, are unphysical and are probably caused by lower signal-to-noise 
near the band edges. 
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Table 1. Targets and Observations 



Source 


Q 

(J2000) 


5 

(J2000) 


d 

(PC) 


Observation Date 
CUT) 


S pro j 

(m) 


■Bpa 
(°) 


Calibrators 


FU Ori 
V1057 Cyg 
V1515 Cyg 


05 45 22.36 
20 58 53.73 
20 23 48.02 


+09 04 12.30 
+44 15 28.54 
+42 12 25.80 


412 
600 
1000 


25 Nov 2010 
15 July 2009 
21 July 2010 


82-84 
78-81 
78-80 


53-55 
46-53 
45-48 


HD 37147, HD 42807 
HD 192985, HD 195050 
HD 194206, HIP 102667 



Note. — S pro j and Bpa are the projected length and position angle of the KI baseline on th e sky during our observa- 
tions. The distance to FU Ori is based on trigonometric parallax measured to the Orion region iReid et al.ll2009h. V1057 
Cyg lies in the North American Nebula, for which the distance has been estimated to be 600 pc jLaugal^ ^Tst raizvsl 
2002). V1515 Cyg resides in the Cyg Rl region, for which the distance has been estimated as 1000 pc llRacineill96*slv 



Table 2. Disk Model Fits 





FU Ori 


V1057 Cyg 


V1515 Cyg 


Xr 


0.5 


1.9 


2.6 


Rm (AU) 


0.07± 0.02 


0.03± 0.07 


0.03± 0.11 


fiout (AU) 


0.49± 0.03 


3.54± 3.14 


1.02± 0.20 


a 


-0.95± 0.10 


-0.50± 0.03 


-0.72± 0.13 


T in (K) 


5529± 779 


4054± 1169 


5973± 1548 


/ 


0.02±0.01 


O.lldbO.01 


0.03±0.03 


A8 (mas) 


0±0 


36± 104 


0± 


Xy2 


l.lrtO.l 


l.OiO.l 


l.Oil.O 


Aniix 


l.OrtO.l 


l.OiO.l 


l.Oil.O 



Note. — Parameters are defined in $H While we also 
fitted inclination and position angle for V1057 Cyg and 
V1515 Cyg, these very poorly constrained by our data 
and so we have not listed the fitted values here. 



